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ABSTRACT 


The  first  section  of  this  paper  uses  the  method  of  stationary 
phase  In  order  to  derive  an  asymptotic  expression  for  the  cross¬ 
correlation  function  for  two  FM  signals  whose  Instantaneous  frequency 
curves  Intersect  at  one  point  In  time-frequency  space.  Assuming  slow 
variation  In  the  amplitude  and  frequency  modulation  functions  allows 
a  simple  geometric  Interpretation  for  the  asymptotic  result  which  Is 
the  square-root  of  an  area  measured  In  time-frequency  space.  This  Is 
shown  In  section  two.  The  case  for  multiple  Intersections  Is  discussed 
and  shown,  through  an  example,  to  depend  on  the  relative  phases  of  the 
Individual  contributions.  This  Intuitive  method  of  evaluating  cor¬ 
relation  functions  explains  why  V-chlrp  and  SQFM  signals  have  pedestals 
in  their  ambiguity  functions. 
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I.  Introduction . 

Frequency  Modulated  (FM>  signals  with  large  time-bandwidth 
products  are  used  in  many  areas,  such  as  Sonar,  Radar  and  Spread 

Oj*) 

Spectrum  Communications.  In  order  to  take  advantage  o-F  the  pulse 

compression  properties  o-f  these  signals,  many  o-f  these 

; 

applications  require  coherent  or  matched  filter  processing  as  the 
appropriate  receiver.  The  receiver  response  is  usually  expressed 
in  terms  of  the  ambiguity  functions  for  the  signals  involved. 
Although  there  are  many  types  of  FM  signals  which  have  pulse 

y 

compression  properties,  th^ir  ambiguity  functions  can  have  quite 
different  properties.  For  example,  linear  FM  signals  (chirp) 

have  a  single  ridge  with  no  pedestal,  whereas  V-chirp  FM  signals 

; 

have  two  small  ridges  near  the  time-frequency  origin  and  an 
additional  pedestal.  Thus, it  is  desireable  to  have  simple  and 
intuitive  methods  for  evaluating  these  ambiguity  functions.  — 

In  the  first  part  of  this  paper  we  use  the  method  of 
stationary  phase  to  derive  an  asymptotic  formula  for  the 

crosscorrel lation  function  of  two  FM  signals  whose  instantaneous 
frequency  curves  cross  at  only  one  point  in  time.  For  this  case, 
a  simple  geometric  procedure  is  provided  for  estimating  the 
crosscorrelation  function,  based  upon  the  overlapping  area  of  two 
templates  generated  in  time— frequency  space.  For  signals  with 
more  then  one  crossing  in  the  time-frequency  plane,  the  geometric 
procedure  provides  an  approximate  upper  bound  on  the 

crosBcorrealtion  function,  since  the  absolute  phase  shift  between 
the  two  signals  may  not  be  known.  In  addition  a  predictable 
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oscillation  in  time  occurs  in  the  crosscorrelation  function  due 
to  the  changing  phase  between  the  components.  Since  the  ambiguity 
function  is  a  crosscorrelation  function,  we  can  apply  these 
results  directly  to  obtain  the  desired  properties. 
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II.  Asymptotic  Theory. 

In  this  section  we  consider  the  problem  of  asymptotically 
estimating  the  crosscorrelation  function  between  two  R1  signal, 
based  upon  the  crossings  of  their  instantaneous  frequency  curves. 
We  shall  deal  with  complex  analytic  signals.  In  the  appendix  we 
establish  the  relationship  between  the  crosscorrelation  function 
of  the  real  parts  of  the  analytic  signal  and  those  for  the 
analytic  signals  themselves. 

We  begin  by  defining  instantaneous  frequency.  Suppose  we  have  an 
analytic  signal  of  the  form 


■3 


By  analytic  we  mean  that  the  Fourier  transform  of  f(t),  F(w), 
satisfies  the  conditicsn 


F(u.)  =0^  «><  o  . 


Thus  the  signal  has  a  one-sided  spectrum.  The  instantaneous 
frequency  for  f(t)  is  the  time  rate-of -change  of  its  phase,  thus 


Ada) 


Hence,  for  example,  if 
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©(t^  =  ^  -fc”  ♦ 


(*) 


for  a  typical  1 inoar  FM  signal,  then 

(0^(0  =  8  *b  **^0  . 


65-) 


We  assume  that  the  two  signals  of  interest  are  of  the  same  form: 


/  =  aw)  e 


—  D  tt)  ft  j 


a-W 


and  the  crosscorrelation  is  given  by 


]■  r>  J* 


(») 


Without  loss  c^F  generality  we  assume  the  analytic  signals  have 
unit  energy.  Thus 


Uiiol*  i-t  =  [  =  1 


Cf) 


By  the  Schwarz  Inequality,  the  maximum  values  of  their  auto  and 
crosscorrelation  functions  are  unity.  Their  instantaneous 
frequencies  are  thus 


a-fr 


and 
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The  two  constants  ,  A  and  B,  are  assumed  to  be  large. 

Substituting  equation  (6)  into  equation  (7)  yields 

(leb 

where 

kl-tjtir  AC4j  k  t-t+ti  ^ 

If  we  assume  that  A  >  B,  without  loss  of  generality,  then  we  have 

.  Qz.) 

Cs) 

Applying  the  method  of  stationary  phase  to  the  integral  in 
equation  (4>  we  have,  as  an  asymptotic  expression 

y.  .  >, 

k ( ^  1^'  3 

where  the  stationary  point,  is  the  unique  solution  of  the 
equation 

W)j  (?)-  ^  • 


It  is  assumed  .hat  g(t)  is  continuous  and  h(t>  is  twice 
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continuously  differentiable  with 


>  o 


Let  us  do  an  example  of  how  this  procedure  can  be  applied. 

We  will  first  consider  the  case  of  the  crosscorrealation  between 
two  linear  FM  signals,  one  an  up  chirp  and  the  other  a  down  chirp. 
We  let 


c  ^ 


o  ±  ^  i  -r 


The  intantaneous  frequency  curves  for  the  two  signals  are 
shown  in  Figure  1. 

The  unit  energy  condition  reveals  that 


f it  «  i-t  =  /  . 


r.tK 


W 

Solving  for  the  stationary  point,  X  =  A,  we  have  that 

A 


and  differentiating  yields 


i-t 


**  A*  ^ 
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A 

Setting  the  derivative  of  h(t>  equal  to  zero  and  solving  for  t,  we 
obtain 


A  T- 1 
=  Z.  * 


The  second  derivative  condition  is  satisfied  since 


At' 


Z  >  o 


J 


^  A 

Solving  for  k<t,T>  and  h(t,T>  we  have 


Thus  we  have  as  an  asymptotic  estimate  of  the  crosscorrelation 


If  we  further  assume  that  the  amplitudes  have  the  form 

UO  =  ^  j  O  ±  -t  i  T-  . 

Thus  the  signals  have  rectangular  envelopes.  Then  in  the 
overlapping  region  <-T  <  '2'  <  T)  we  have 
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'A 


We  observe  that  the  total  frequency  deviation  is  AT 
radians/second  or  radians/second,  then  Me  obtain  that 

^  *T*  \  ^  J  *"  Ta.F^ 

We  see  that,  as  the  time-bandMidth  product  gets  large,  the 
crosscorrelation  function  diminishes  as  the  reciprocal  square- 
root  of  the  time-bandwidth  product. 

In  order  to  check  the  asymptotic  results,  several  computer 
examples  were  run.  The  signals  were  of  the  form 

Z 

F(I)  =  A(I)  EXP(BI  *  WI>,  I  s 

The  ampitudes  were  rectangular  except  for  a  lOX  raised  cosine  at 
the  ends  of  the  signals.  Figures  2.  and  3.  show  the  two  signals 
for  N  —  128.  The  real  and  imaginary  parts  of  the  complex  signals 
are  plotted  consecutively  on  each  graph.  The  parameter  are 

a)  Upsweep  B  »  0.003,  W  »  o.3 

b)  Downsweep  B  =  -0.003,  W  =  1.068 

Calculating  the  time  duration  and  bandwidth  we  have 

Time  duration  »  (0.9) C128)  «  115.2 

Bandwidth  »  (0.9X0.768/2  )  -  0.110 


y.f  y>y. 
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Figure  2,  Real  and  Imaginary  parts  of  128  point 
Linear  FM  Upsweep  signal. 
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Thus  the  crosscorrelation  estieete  is 


«  0.19S6 


Figure  4.  shoMS  the  coeputed  crosscorrelation  function  overlayed 
with  the  asymptotic  estimate,  shown  as  a  staight  line. 

Suppose  the  lOX  raised  cosine  window  is  not  included  then 

Time  duration  »  128 

Bandwidth  »  0.768/2 

and  the  crosscorrelation  estimate  is  0.179.  The  results  of  this 
case  are  shown  in  Figure  S. ,  where  the  computed  crosscorrelation 
function  and  the  estimate  or  overlayed. 

A  seconds  example  with  larger  time -bandwidth  product  is  given 
below.  Here  the  parameters  are 


N  «  2S6 

a)  Upsweep  B  »  0.003,  W  »  0.3 

b)  Downsweep  B  a  -0.003,  W  «  1.836 

Figure  6.  shows  the  first  half  of  the  crosscorrelation  functions 
of  the  signals  overlayed  with  the  asymptotic  estimates.  The 
results  shown  are  for  signals  with  and  without  the  10%  raised 


cosine  window 
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III.  Geometric  Interpretation. 

In  this  section  we  establish  that  a  simple  and  intuitive 
geometric  interpretation  can  be  applied  to  the  asymptotic  results 
for  the  magnitude  of  the  crosscorrelation  function  in  equation 
(13).  In  particular  we  show  that  the  magnitude  of  the  right  side 
of  equation  (13),  the  asymptotic  estimate  of  the  crosscorrelation 
function,  is  obtained  directly  from  the  area  of  a  parallelogram 
in  the  time-frequency  plane. 

We  begin  by  defining  a  par al 1  el ogr am  in  terms  of  two  sets  of 
parallel  lines  in  the  time-frequency  plane.  Figure  7.  shows  the  4 
lines.  The  area  of  the  parallelogram  is 


Area  =  xy  sin(^). 


O' S') 


We  must  calculate  the  three  quantities  separately.  To  obtain 
si n (  )  we  observe  that 

tan  =  A,  tan  y*  =  B. 

Note  that  w  is  defined  with  a  negative  slope.  Thus  we  have  that 

tan(90  -«C>  =  1/A,  tan(90  -/ff  >  =  1/B. 


Taking  inverse  tangents  and  adding  we  have 


tan*'  (1/A)  +  tan**  (1/B)  =  ISO  -  -  P  =  6 


Figure  7.  Geometry  for  Calculating  the  Area  of  Parallelogram 
Enclosed  by  Two  Sets  of  Parallel  Lines. 
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Thus 


tan  d  —  tan  1  tan  <1/A> 


"a6  -r 


tan  (1/B> 


r- 


and  from  this  Me  obtain 


sin  9  s 


A  -t  B 


^(AViU  bS. 


To  obtain  x  and  y  we  recall  that  the  diameter  of  the 
inscribed  circle  for  a  triangle  is  calculated  as  the  ratio  of  any 
side  to  the  sin  of  the  opposite  angle.  Thus 

c  .  X _  jD_  «  y  . 

SI  n  B  “  s  in  Pk  ^  Sin  ©  Sia  /J 


C  S  in  < 


D  Sin  p 

S  in  © 


O'?) 


s  lA  = 


Thus  combining  (16),  (17),  (18)  and  substituting  into  equation 

(IS),  we  obtain  that 
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If  we  now  assume  that  the  the  two  amplitude  functions,  a(1:)  and 
b(t),  are  slowly  varying,  and  that  the  two  instantaneous  frequency 
functions  are  slowly  varying,  then  we  can  associate  the  quantities 
in  equation  (20)  with  those  in  the  asymptotic  equation  (22).  Here 
the  amplitudes  are  assumed  to  be  constant,  whereas  the  the 

frequency  functions,  w^  and  w^,  are  replaced  by  the  first  two 

^  A 

terms  of  their  Taylor  series  expansions  about  the  point  t,  for 

fixed  T  •  Then  we  may  observe  that 

a.  A 

Ds  b 

j 


and 


8=  - 


,  i  «>3  Lt, 
ztK  i  t 


Then  we  have  finally  that 


\  <o|  ~ 


'/*. 


and  thus  the  crosscorrelation  function  can  be  obtained  simply  be 
taking  the  square— root  of  the  parallelogram  described  above. 

The  results  are  true  for  the  two  linearly  frequency  modulated  signa! 
provided  in  the  example,  since  all  the  assumptions,  requiring  slow 
variations,  are  true.  Although  we  have  not  established  these 
results  when  the  slow  variation  assumptions  are  loosened.  Several 
computer  examples  indicate  that  the  results  are  still  approx i mat ly 
true,  in  those  regions  of  X  which  the  assumptions  are  not 
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true. 

We  can  now  see  that,  in  order  to  obtain  the  asymptotic 
estimate  of  the  crosscorrelation  function,  one  need  only  construct 
a  template  for  each  signal,  such  as  seen  in  Figure  ,  and  overlay 
the  two  templates  for  a  fixed  value  of  ,  and  take  the  square-root 
of  the  reulting  common  area.  Moving  one  template  horizontally  with 
respect  to  the  other  would  yield  a  new  value  of  crosscf3rrelation. 
Moving  the  same  template  vertically,  simulating  a  narrowband 
Doppler  shift,  would  provide  a  value  of  the  narrowband  cross¬ 
ambiguity  function  for  the  two  signals.  The  time  and  frequency 
coordinates  at  which  we  are  evaluating  the  crossambigui ty  function 
is  the  horizontal  and  vertical  shifts  respectively. 

Although  it  may  seem  confusing  as  to  why  we  plot  the 

*  4 

amplitude-squared  along  the  frequency  axis,  i.^.  L*  »  a  Ct,f>, 

observe  that  for  the  rectangulare  amplitude  case,  as  in  the 

«  K  ^ 

examples  shown  in  the  previous  section,  a  (t,T>  =  1/T,  which  has 
the  units  of  frequency,  so  that  the  reults  are,  at  least. 


consistent. 
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IV.  Multiple  Intersections. 

The  case  of  multipe  intersections  of  the  instantaneous 
frequency  curves  for  the  two  signals  is  considerably  more 
complicated.  The  reason  for  this  complexity  becomes  clear  when 
we  consider  a  simplified  case.  Suppose  that  for  a  particular 
delay,  the  instantaneous  frequency  curves  have  two  intersections. 
Then  we  may  calculate  the  crosscorrelation  contribution  from  each 
intersection.  Now  the  crosscorrelation  function  has  two 
stationary  points  and,  hence,  has  two  complex  terms  which  must  be 
summed  to  obtain  the  complete  result.  In  principle  we  have  the 
asymptotic  expression  for  the  phase  in  equation  (13)  and  could, 
thus,  calculate  the  complex  sum  of  the  two  contributions.  This 
would  only  be  possible  if  the  two  signals  were  specified 
analytically.  However,  if  the  signals  were  specified  by  their 
instantaneous  frequency  plots  and  amplitudes,  then  we  do  not  have 
any  knowledge  of  the  absolute  phase  relationship  between  the  two 
contributions.  Since  the  addition  of  the  contributions  is 
complex,  the  resulting  magnitude  could  have  any  value  between  the 
sum  and  difference  of  the  individual  contributions,  depending  upon 
their  phase  relationship. 

An  even  further  complication  is  that,  the  phase  of  each 
contibution  is  time-dependant.  Thus  the  complex  addition  can 
change  for  each  crosscorrelation  estimate.  This  would  be 
evidenced  by  an  oscillation  in  the  crosscorrelation  function 
between  the  largest  and  smallest  values  of  the  complex  addition. 

In  order  to  demonstrate  this  effect,  Figure  9.  shows  the 
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crosscorrelation  -function  for  two  signal  with  several  intersecting 
regions.  The  signals  are  128  points  in  duration  with  maxi mum 
frequency  of  2^  .  Both  signal  have  10%  raised  cosine  windows  on 
each  end.  The  first  signal  is  the  sum  of  two  linear  Fli  signals 
with  frequency  ranges  <. 5,1.5)  and  (2,3),  while  the  second 
signal  is  a  single  linear  Fii  sweep  with  range  (3,1.25).  If  one 
constructs  the  two  templates  and  passes  them  over  one  another 
horizontally  then  we  observe  that  first  there  is  no  overlap,  then 
a  small  region  of  one  intersection,  a  region  of  two  intersections, 
and  finally  a  region  of  one  intersection.  The  region  of  two 
intersections  shows  a  rapid  oscillation  due  to  the  relative 
phase,  which  is  changing  in  time. 
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COBR  Crosscorr'e I a-t i on  Function  Example 


Figure  9.  Crosscorrelation  Function  Example  with 
Several  Intersecting  Regions. 


V.  Conclusions 


We  have  shown  that  a  simple  and  intuitive  geometric 
procedure  can  be  used  to  obtain  an  asymptotic  estimate  of 
crosscorrelation  functions  for  FM  signals  whose  instantaneous 
frequency  curves  cross  at  only  one  point.  In  other  casc»s  this 
procedure  provides  an  approximate  upper  bound  on  the  behavior  of 
the  crosscorrelation  function,  in  regions  of  more  than  one 
intersecting  region.  The  procedure  can  be  applied  directly  to 
narrowband  ambiguity  functions  since  the  a  Doppler  shift  is 
simpley  a  vertical  shifting  of  the  overlapping  template.  We  can 
easily  see  reason  why  the  autoambiguity  function  for  a  Vchirp  or 
SQFM  signal  has  a  pronouced  pedestal.  Both  of  these  signal  pass 
through  the  same  range  of  frequencies  twice.  Hence,  for  any 
shift,  there  will  always  be  one  intersection  of  the  instantaneous 
frequency  curves.  For  the  Vchirp  signal  the  pedestal  is  completly 
flat,  since  the  parallelogram  is  the  same  area  for  all  shifts. 
In  the  case  of  the  SQFM  signal  the  parallelogram  reduces  area  as 
you  move  from  the  origin. 

Perhaps  the  most  useful  aspect  of  this  geometric  procedure 
is  the  intuitive  ability  for  quickly  evaluating  the  properties  of 
crosscorrelation  and  ambiguity  functions  of  FM  signals.  Further, 
if  one  must  evaluate  the  properties  of  many  FM  signals,  a  template 
need  only  be  constructed  for  each  one.  Thus,  one  can  evaluate 
these  functions  rather  quickly. 

Finally,  as  new  electronic  signal  proccessing  devices,  which 
can  represent  signals  in  time-frequency  space,  become  available. 
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Appendix. 

In  this  appendix  we  establish  that  if  we  have  two  analytic 
signals,  k(t)  and  h(t),  then  the  crosscorrelation  of  their  real 
parts  is  equal  to  the  real  part  of  their  crosscorrelation 
function.  Further  that  the  magnitude  of  the  complex 
crosscorrelation  function  is  the  envelope  of  the  real 
cr osscor r el at i on . 

To  begin,  we  assume  that  f ^ (t>  and  f^(t>  are  real,  unit 
energy  signals,  and  the  transform 


A 

r  1 

1  ( 

f  <t)  = 

H  [f  (t)| 

=  -—  If  <r>dr /(t- 
) 

^  ao 

t ) 

(A-1) 

the  Hilbert 

transform. 

Then  we  form  the 

analytic 

signals 

^  k<t) 

h<t) 

to 

Denoting  the  inner  product  as 

<k,h)  =  j  k<t)h^<t)dt  (A-3) 

and  using  the  properties  of  Hilbert  transforms  and  their  spectra 
it  can  be  shown  that  if 

F,  <w)  =  A.(w)  e  (A-4-a) 


=  J^f,<t)  +  <A“2-a> 

-  [  f»<t) 
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Fj^(w}  =  Aj|.(w)  a 


(A~4-b) 


then 


(k,h> 


=  <f,  .f,>  ♦  I  A,<->ftj 


(w>Bin  (aCCw)  )dw  /  (A— 5) 


•] 


where 


0^(w)  =  0|  (w)  -  (w) . 


Recognizing  that  if 


f  I  <t) 

f^<t) 


f  <t) 

g<t  +t) 


Then  the  inner  product  is  a  crosscorrelation  function  and 


Further  since  has  a  one  sided  spectrum  then  its  real  and 


imaginary  parts  are  themselves  Hilbert  transforms  and  hence 
iBbuCOl  envelope  of 
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